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Abstract 

The complete two-loop expression for the jet function J{p^,fi) of soft-colUnear effective theory 
is presented, including non-logarithmic terms. Combined with our previous calculation of the 
soft function S{a>,iu), this result provides the basis for a calculation of the effect of a photon- 
energy cut in the measurement of the B X^y decay rate at next-to-next-to-leading order in 
renormalization-group improved perturbation theory. The jet function is also relevant to the 
resummation of Sudakov logarithms in other hard QCD processes. 



1 Introduction 



A significant effort is currently underway to complete the Standard Model calculation of the B 
X^y decay rate at next-to-next-to-leading order (NNLO) in renormalization-group improved pertur- 
bation theory. It is motivated by the fact that the relatively large branching ratio for this decay, 
combined with the increased precision in its measurement at the B factories, make this an excellent 
way to probe for hints of new flavor physics. The experimental detection of B ^ X^y events relies 
on the reconstruction of a high-energy photon, whose energy in the 5-meson rest frame exceeds a 
value £0 ~ 1.8 GeV. The theoretical analysis of the partial inclusive B X^y decay rate with a cut 
Ey > Eq must deal with short-distance contributions associated with three different mass scales: the 
hard scale m/„ an intermediate scale ^fri%K, and a soft scale A, where A = m/, - 2Eq ^ 1 GeV Q. 
The cut-dependent effects are described in terms of two perturbative objects called the jet function 
and the soft function, which for an analysis of the decay rate at NNLO are required with two-loop 
accuracy. The two-loop calculation of the soft function has been presented in 1 2 1 , while that of the jet 
function is described in the present work. As a by-product of our analysis we calculate the two-loop 
anomalous-dimension kernel of the jet function. 

The jet function j(L,/j.) needed in the factorization formula for the partial B Xgy decay rate 
is related to the original jet function J{p^,/u) of a massless quark in QCD [|4i| by 

j(\n^,iu)^ dp^J(p\lu). (1) 
^ /^^ ^ Jo 

While the perturbative expression for J{p^,/j.) involves singular distributions (see, e.g., lElllIl), the 
function j has a double-logarithmic expansion of the form 

;(L,^) = 1 + 1; (^)" (^f > + b^L + • • • + b^l.L^"-^ + Z^^jL-) . (2) 

n=l ^ ^ 

By solving the renormalization-group equation for the jet function order by order in perturbation 
theory, the coefficients Zj^'^q of the logarithmic terms in @ can be obtained from the expansion co- 
efficients of the jet-function anomalous dimension and the /3-function, together with the coefficients 

arising in lower orders [3|. The two-loop calculation performed in the present work gives the 
constant b''^^ and provides the first direct calculation of the two-loop anomalous dimension of the jet 
function. We also note that from our result for j(L,/j.) one can derive the two-loop expression for 
J(p^,fj.) in terms of so-called star distributions 

We stress that even though our primary goal is to improve the theoretical analysis of 5 — > X^y 
decay, applications of our results are not confined to flavor physics. Indeed, the jet function is 
a universal object, which enters in many applications of perturbative QCD to jet physics, deep- 
inelastic scattering, and other hard processes. The two-loop calculation of the function j{L,^) is 
described in Sectional It follows closely our calculation of the soft function in however, the 
evaluation of the two-loop master integrals is considerably more complicated in the present case. In 
Section|3]we briefly discuss jet-function moments and their renormalization-group evolution. 
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2 Two-loop calculation of the jet function 



The factorization properties of decay rates and cross sections for processes involving hard, soft, 
and collinear degrees of freedom become most transparent if an effective field theory is employed 
to disentangle the contributions associated with these different momentum regions. Soft-coUinear 
effective theory (SCET) has been designed to accomplish this task ||Hl|51[TnilIIl. In the context of 
SCET the jet function is defined in terms of the hard-coUinear quark propagator |l6l|5| 

'^n-p Jip\fi) = J d'x e-'"-' <0 1 T [XhM %c{x)] 1 0) , (3) 

where /i is the renormalization scale, and n and n are two light-like vectors satisfying n ■ n = 2. 
For simplicity we suppress color indices on the quark fields. The propagator is proportional to a 
unit matrix in color space. The composite field Xf,c(x) = Sl(x-)Wl^(x)^(x) ifTTl [12, 131 is the 
gauge-invariant (under both soft and hard-coUinear gauge transformations) effective-theory field for 
a massless quark after a decoupling transformation has been applied, which removes the interactions 
of soft gluons with hard-coUinear fields in the leading-order SCET Lagrangian ['^l. In the absence of 
such interactions the hard-coUinear Lagrangian is equivalent to the conventional QCD Lagrangian, 
and we can rewrite the propagator in terms of standard QCD fields as 

'^n-pj(p\iu) = J J^xe-'>^<OIT|^W+(0)«A(0)^WWW^||0). (4) 

The quark fields are multiplied by Wilson lines 

W(x) = P exp [ig J dsn ■ A{x + sn)j , (5) 

which render the expression © gauge invariant. Note that the Wilson lines are absent in the light- 
cone gauge n-A = 0. For this reason the function ^ is sometimes referred to as the quark propagator 
in axial gauge. Lorentz invariance dictates that the QCD propagator in the presence of these Wilson 
lines contains two Dirac structures proportional to / and $. The Dirac matrices appearing to the left 
and right of the field operators in (|H) project out the terms proportional to The jet function J is 
the discontinuity of the propagator, i.e. 

Jip\n) = ^ Im [iJip\/^)] = 5{p^) + 0{as) . (6) 
Finally, we calculate the function j from the contour integral 

j(ln^,//) = dp^J(p\M) = f dp^J{p\p). (7) 

Ip2|=e2 

Our calculation of the jet function employs the representation Q of the function J^{p^,iJ) in 
terms of ordinary QCD quark and gluon fields. The relevant two-loop diagrams are shown in Fig- 
ureHl Equally well, one could use the SCET Lagrangian together with © to perform the calculation. 
In this case diagrams in which a quark emits more than one gluon at the same vertex would also be 
present, in addition to the topologies shown in FigureUl Also, the analysis would be complicated by 
the fact that the SCET Feynman rules are more complicated that those of QCD. 
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Figure 1: Two-loop diagrams contributing to the jet function in QCD. Gluons emitted from the 
crossed circles originate from the Wilson lines. Not shown are additional diagrams resulting from 
mirror images in which the two external points are exchanged. The first diagram is the full fermion 
two-point function, not just the one-particle irreducible part. 

2.1 Evaluation of the two-loop diagrams 

We first discuss the evaluation of the bare quantity jbare(2^) and later perform its renormalization. 
Let us begin by quoting the result for the one-loop master integral 



d'^k = ini i-p^ - (n ■ pT' J{a, b, c) , (8) 

with 

- b)T(^ - a - c)Y(a + b - ^) 
J(a, b, c) = — ^ . (9) 

At two-loop order, the most general integral we need is (omitting the "-l-zO" terms for brevity) 



-a~h-c 



(_ \ \-ai-a2-ai-bi~h2~b}-ci-C2 



{Bf ipf' [{k - lyf [{k + pft [(/ + pff' [{k + l + pff' (n ■ kr in ■ If' 

1/ r.\d-a\-a2-ai-b\-b2-h'i _ _ _ 

= -n"\-p^) (n-p) ' Jiai,a2,a3,bi,b2,b3,Ci,C2). (10) 



We use the same standard reduction techniques as in the two-loop calculation of the soft function Q 
to express all integrals we need for the evaluation of the diagrams in Figure[l]in terms of four master 
integrals M„. Introducing the dimensional regulator e = 2 - d/2, we obtain 

Ml = 7(1 , 1 , 0, 0, 0, 1 , 0, 0) = ^ ~ f ^ 

r(3 - 36) 
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M2 = 7(1, 1,0, 1,1, 0,0,0) = 7(1, 1,0)2 = 
M3 = 7(1,0,1,1,1,0,0,1) 



r(2 - 26) 



r(26) r\ r\ {xxy^n 



Jdx \ dy^ 
Jo 1 - xy 

M, = 7(1, 1,0, 1, 1, 1, 1, 1) = e-^'^^ ~ T ^ ^ ^ ' ^^^^ 

where we use the shorthand notation x = \ - x and y = 1 - The evaluation of the first three 
integrals is straightforward. In the case of Ma the double parameter integral resulting from the loop 
integrations can be expanded in e without difficulty. The calculation of the master integral A/4, which 
is needed for the evaluation of the seventh graph in Figure [H is notably more complicated. 

To tackle this last integral we use the Mellin-Bames technique ifT^lfTSl . The basic strategy is to 
first introduce Feynman parameters to perform the loop integration and then introduce Mellin-Bames 
parameters to carry out the Feynman parameter integrations. What makes this method powerful is 
that (after analytic continuation to 6 « 0) the Mellin-Bames integrands can be Taylor expanded 
about 6 = 0. To start, note that after performing the loop integral over k using conventional Feynman 
parameters the result for Ma, can be written as 



M4 = in-'^ (-p^f''' (n ■ pf r(3 - f ) J d'^l 



1 



f {l + pfn-l 



X 



{ dx { dy i ^- J. (12) 

Jo Jo n-p+yn-l / ^-2^72 _ , „y _ yjcMn^^"^ 



(-x^yyP - xxy(l + pY - xxyp^) 
We now introduce two Mellin-Bames parameters via 



1 r Ti-w) Via + w) 

27r« Jc-ioo r(a) 

to break apart the last denominator in dT^ into a product with factors f-, (/ -1- p)^, and p^. (We do 
not introduce a Mellin-Barnes parameter for the denominator {n ■ p + y il ■ I), as this would lead to 
ambiguities in the treatment of poles in the resulting light-cone propagators.) We then introduce ad- 
ditional Feynman parameters and perform the loop integration over I. We use a third Mellin-Barnes 
parameter to simplify the resulting Feynman parameter integrals, which can then be expressed in 
terms of F-f unctions. This leaves us with the following three-fold Mellin-Bames integral: 

M4 = — , Tl r^d.^ r(-w,)r(-.,)r(-.3)r(i ..0 



(ini)^ \j j,,_,.eo r(-26) r(i - 26 + wi + W2 + W3) r(i - W2) m - w^) 
X r(-6 - W2) r(6 -W2- W3) r(i - 6 -1- wi -1- W3) r(i + e + W2 + W3) r(i + wi + W2 + w^) 
X r(w2 - 6)r(-i - 6 - wi - W3) . (14) 
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In deriving this representation we have interchanged loop, Feynman, and Mellin-Barnes integrations. 
Careful inspection reveals that the representation (fT?l) is valid if the real parts of all F-functions are 
positive, i.e., if the poles of each F-function are either all to the left or all to the right of the integration 
contours. With the choice ci = -|,C2 = -|, and C3 = -|, this condition is fulfilled if-|<6<-|. 
Starting in the allowed range and increasing e, we see that at 6 = -| the first pole in F(- 1 -e-wi -W3) 
crosses the contour from the right to the left. At e = the first pole in F(w2 - e) crosses from the 
left to the right. The arguments of all other F-functions remain positive up to 6 < |. To obtain a 
representation that is valid around e = 0, one has to either deform the contours such that the crossings 
are avoided, or separately take into account the contributions of the poles that end up on the wrong 
side of the contours as 6 ^ 0, as proposed in 1 15 1. The residues of these poles have again the form 
of Mellin-Barnes integrals, however with one integration less than the original expression (fT4l) . One 
analyzes these contributions with the same method as the original integral and continues until one 
ends up with a representation which is valid around e = 0. Once this is achieved, the integrands are 
expanded in e, the integration contours are closed, and the integrations are rewritten as sums over 
the residues of the poles. In fact, since the original integral includes a factor 1 /Y{-2e) = 0{e), only 
contributions which arise from poles that cross a contour need to be included in the limit 6^0. 
Very recently, the continuation in e with subsequent numerical evaluation of the integrals has been 
automatized lITClfTTll . We have used the public code of [il7il to check our analytical result for M4. 

As a further independent check of our result we have numerically evaluated the multi-dimen- 
sional integral over Feynman parameters, which is obtained after performing the loop integration 
over / in dT^ using conventional methods. As is evident from dTTT) the integral M4 is divergent, and 
the divergences need to be isolated in order to perform the numerical evaluation. We use the method 
of sector decomposition ifTSl fT9l l20ll . which allows one to systematically disentangle overlapping 
singularities in Feynman integrals. This procedure splits the integral into a large number of terms 
in which all singularities are factorized. Because it leads to large algebraic expressions, the sector 
decomposition is performed using computer algebra. After numerical integration of the resulting 
expressions we reproduce the analytical result for the integral M4 with a numerical precision of 
better than 1 part in 10^. 

With the integrals at hand, the evaluation of the two-loop diagrams in Figure[T]is straightforward. 
We write each diagram as a sum of integrals of the form (fTOI) and express those in terms of the four 
master integrals Mi, ... , M4. Summing up the results for the individual diagrams, we obtain 



7bare(e') = l + ^l^l O 



An \ff 



4 3 _ 2 /- . 37r2 28 



_,_,V-.^.^14----^. 



+ (^) (^j Cf [CpRpie) + CAKAie) + TpnfKf{e)\ + . . . , (15) 

Z„ = 1 -ySo — + • • ■ = 1 - I — Ca - - TpTlA — + .... (16) 

Ane \ 3 3 ^/ Ane ^ ' 

Here = asip.) is the renormalized coupling. Note that the bare jet function is scale indepen- 
dent, since the bare coupling Z^as/i^'^ does not depend on the renormalization scale. The two-loop 



where 
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coefficients are 



... X 8 12 (65 16n^\ 1 /311 ^ , 124 \ 1 

1437 _ 301^ ^n3^_ 

8 12 90 V 
, , 11 /233 1 /4541 SS/r^ \ ^ 

86393 1129;r2 IVtt^ 514 

+ ^"3 + 0(6) , 

648 108 180 9 ^ 
4 38 / 373 107r2\ 1 7081 957^^ 128 

We have checked that the divergences of the first diagram in Figured] (together with the appropriate 
one-loop counter term, and accounting for the renormalization of the gauge parameter) reproduce 
the known result for the two-loop quark wave-function renormalization ED • A stringent check of 
the remaining diagrams with gluon emissions from the Wilson lines is that their divergences must 
cancel against the jet-function renormalization factor, which we will now derive. 

2.2 Renormalization of the jet function 

The renormalization of the bare jet function proceeds in complete analogy to that of the bare soft 
function discussed in Q, to which we refer the reader for a more detailed discussion. The procedure 
is more complicated than in conventional applications of renormalization owing to the presence of 
Sudakov double logarithms. 

We define an operator renormalization factor for the jet function via 

2 .,\ / j>2 7r„1 „'2 ,A T 111 



J{p\^i) = J dp"Z{p\p'\fi)AUp'l, (18) 

where Z absorbs the UV divergences of the bare jet function, such that the renormalized jet function 
is finite in the limit 6 ^ 0. In the MS regularization scheme, we have 

Zip\p'\iu) = 6ip^-p'^) + J] -Z^'\p\p'\n). (19) 

k=\ ^ 

The relations 

rjet = 2a, -— , 2a, — = 2a, —— ® Z^ + fiia,) —— + —— , (20) 

oa, oa, oa, oa, olnyu 

with n > \, connect the coefficient of the 1/e pole in the Z factor to the anomalous dimension and 
moreover imply a set of consistency conditions among the coefficients of the higher pole terms. The 
symbol ® represents a convolution in p^. The term dZ^"^^ IdXnjj. arises because the Z factor depends 
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both implicitly (via the renormalized coupling constant) and explicitly (via Sudakov logarithms con- 
tained in star distributions) on the renormalization scale 0- 

To all orders in perturbation theory the anomalous-dimension kernel of the jet function has the 
form 



rjet(p^p'^)") = 2rcusp(Q',) 



1 



p2 _ pfl 



+ 2y(a,,)6ip'-p'^). 



(21) 



where Fcusp is the cusp anomalous dimension associated with the Sudakov double logarithms, while 
y' controls the single-logarithmic evolution of the jet function. The definition of the star distribution 
can be found in The corresponding integro-differential evolution equation reads 



dJ(p^,H) 
Jin// 



cusp 



dp" 
Jo 



2 J(p",iu)-J(p\m) 



p2 _ pr2 



(22) 



We have derived relation (jSTTl by requiring that the B X^y decay rate be renormalization-group 
invariant and using the known evolution equations for the soft function ll22ll and for the hard match- 
ing coefficient IHllHll- Denoting by Z[„] the coefficient of (aJ4n)" in Z{p^, p'^,/u), we obtain from 
(ESI 



Z[0] = <5(P -/'), 
Z, 



-[1] 



-[2] 



6(p^-p'^) 



2e4 



-I- 



+ 



_r| Toiyj - jfio) n 

e^^ 62 ^ 26 



4 12 °) 62 26 



p2 _ p/2 



r2 



In^ 



p2 _ p/2 



(23) 



where the expansion coefficients of the anomalous dimensions and /3- function are defined as 

DO ,1 OO , , 

Z(as\ I lo:s\ 



n=0 

da 



n=0 



' 11=0 



(24) 



The expression for the two-loop cusp anomalous dimension can be found, e.g., in our previous 
paper 1 2 1 . The two-loop anomalous dimension of the jet function has never been calculated directly, 
but it was inferred in 1 1 1 from existing two-loop results for jet-function moments in deep-inelastic 
scattering [.24] . In that way one obtains 



7o = -3Cf , 



(25) 



1769 Il7r2 \ ^ ^ /242 4n^ 

+ 40Z'3 + CpTFUfl + 

54 9 ^7 ^ ^ ^ \ 27 9 



It is straightforward to show that the function j(\n^,fi) obeys the same evolution equation as 
the original jet function J(p^,ij), i.e. 



(26) 



where jbaie(2^) is the quantity we have calculated in Section ITTI Expanding this relation in pertur- 
bation theory we obtain j[o] = j^^^ = 1 and 



i[l] - fu{ + ^[1] 



J[2] - J[2] 



;bare . v 



■bare 



[2] ® J[0] 



(27) 



The first term on the right-hand side in each equation corresponds to the contribution (fT31) obtained 
from the loop diagrams. The remaining terms correspond to the counter- term contributions. Explic- 
itly, we find 



•C.T. 

^[1] 



•C.T. 

J [2] 



— r + — + — In — , 



1^ + 1^ + liillnC 



-I- 



6 6 



+ 



je - (A(l - e) 



n- 



12 



1 y\ 
— + 



+ 



n , ro(yo^ - 

63 62 



vSo) 



2e 



i„4 



(28) 



Together with the results for the bare one-loop jet function from ( [T51 this yields explicit expressions 
for the counter terms. When adding these contributions to the bare jet function we find that all 
lie" pole terms cancel, so that the limit 6^0 can now be taken. This implies, in particular, that 
we confirm by direct calculation the expression for the anomalous-dimension coefficient y[ given 
India). 



2.3 Results 

The logarithmic terms in the renormalized jet function have been determined in |T| by solving the 
renormalization-group equation perturbatively. At two-loop order, it was found that 



KUn) = 1 + 



-I- 



An 
An 



(29) 
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Our results for the logarithmic terms agree with the above expression. The one-loop coefficient b)^ 

^0 



(1) 



was derived in 115101. The main new result is h^^ , the constant term at two-loop order. We obtain 



h^^ = {l-n')Cp, 

, (2) _ ^2 /205 ein^ Un' 1 8/ U T T / ^^^^^ ^^^^^ / \ 

I 4057 687r2 16 \ 

It is interesting to compare the exact answer for the coefficient b^^^ with the approximation obtained 
by keeping only the terms of order Pqo]. In the absence of exact two-loop results it is sometimes 
argued that the Poo] terms constitute the dominant part of the complete two-loop correction. In the 
present case, we obtain for Nc = 3 colors ~ -16.25jc - 128.78 ~ -145.04, where x = ^fio = I 
for = 4 light flavors. Keeping only the /Soa^ terms would give -16.25, which is off by an order of 
magnitude. This illustrates the importance of performing exact two-loop calculations. 

We now briefly discuss the impact of our results for phenomenology. Besides the jet function 
j{L,ij) itself, it is useful to consider a related function j{L,iJ.) obtained by replacing the n-th power 
of L in © with an n-th order polynomial, L" In(L), where at two-loop order we need 

hix) = X , hix) = X + —X-2^2,, 

hix) = x^ + ^, Uix) = / + nh^ - 8^3 ^ + ^ ■ (31) 

The function j enters in the factorization formula for the partial B Xsj decay rate with a cut on 
photon energy f3l . For the case of A^c = 3 colors and n/ = 4 light quark flavors we get 

j{L,i^) ^ 1 +(-0.304- 0.3 18L + 0.2 12L2)a,(ju) 

+ (-0.918 + 0.926L + 0.079L2 - 0.114L^ + 0.023L^) a^Ou) + . . . , 
7iL,iu) « 1 + (0.045 - 0.3 18L + 0.212L2) a,(ju) 

+ (-0.185 + 0.145L + 0.301L2 - 0.1 14L^ + 0.023L4) a^,(M) + ... . (32) 

The two-loop corrections to j are very large and, for realistic parameter values, can even dominate 
over the one-loop corrections. However, the two-loop corrections are much smaller for the function 
j. Figure El shows the dependence of the two jet functions on L = ln(2^/;U^) for a fixed scale 
ju ~ 2GeV chosen such that a sip.) = 0.3, corresponding to a renormalization point appropriate for 
the calculation of the partial B X^y decay rate with a cut Ey > 1.8 GeV. The two-loop effects 
calculated in this work impact the jet function j at the 10% level, while their effect on the function 
j is at the level of 2% or less. The latter finding suggests a good convergence of the perturbative 
expansion at the intermediate scale in the analysis of 5 ^ Xsj decay. In addition to the one- and 
two-loop predictions, the figure also displays the results obtained if only terms of order /JqQ;^ are kept 
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Figure 2: One- and two-loop predictions for the jet functions jiL,fi) and j(L,fj.) evaluated at asQu) = 
0.3. The dashed lines show the one-loop results, while the solid lines give the complete two-loop 
results derived in the present work. The gray lines are obtained if only the /3oa^. terms are kept in the 
two-loop contributions. 

in the two-loop coefficients. In both cases this provides a poor approximation to the exact two-loop 
results. We also note that the jet functions by themselves are not renormalization-group invariant, 
so it is meaningless to study their dependence on the scale fi for fixed Q^. In physical results such 
as the expression for the B X^y decay rate and photon-energy moments given in yj|3J, the scale 
dependence of the jet function cancels against that of other renormalization-group functions. 



3 Moments of the jet function 

In the analysis of hard QCD processes such as deep-inelastic scattering it is often convenient to 
introduce moments of the jet function defined as (see e.g. [25 1) 



Jo 



Q 



2 \N-1 

2 



JNiQ\fi)= I dp'\l-^\ Jip^iu). (33) 



Whereas in inclusive B decays the scale setting the upper integration limit in dU) is an intermediate 
(hard-coUinear) scale, «; m^, which is of order the invariant mass squared of the final- state 
hadronic jet, the variable in is set by a characteristic hard scale of the process. In the large-A'^ 
limit the integral receives leading contributions only from the region ~ Q^/N <^ Q^. The scale 
Q^/N is the analog of the intermediate scale in B ^ X^y decay. 

Using an integration by parts, it is straightforward to express the jet-function moments in terms 
of integrals over the function j calculated at two-loop order in the present work. We obtain 

JiiQ\lu) = j(ln^,n), MQ\n) = (N-\) f dx(l-xf~^j(\n^,fi), (34) 

where the second relation holds for N >2.lt follows from (E^ that the moments obey the evolution 
equation 



dlnfj. 



- 2reusp|ln^-//;v-ij + 27^ 
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2Y 



cusp 



Jo 



N-l 



(35) 



where //a^-i = Ylt=\ \ is the harmonic number. These resuks simplify greatly in the large-A^ limit. 
We find that the moments are given by 



(36) 



and that they obey the local evolution equation 



d In lu 



Jj,{Q\^i) + 0\-\. 



1 



(37) 



In deriving this result we have used that the second line in (l35t is suppressed in the large-A^ limit, 
since p^lQ^ = 0(1 /N) in the argument of the logarithm. This local evolution equation can be 
integrated using standard techniques. 



4 Conclusions 

We have calculated the two-loop expression for the jet function j{L,/u) defined in terms of an in- 
tegral over the hard-coUinear quark propagator in soft-coUinear effective theory. This quantity is a 
necessary ingredient for the NNLO evaluation of the B X^y decay rate with a cut on the pho- 
ton energy. Moreover, since the jet function is universal, it appears in many other applications of 
perturbative QCD. The results obtained in the present work, when combined with |2|, provide a 
complete description of low-scale effects in the analysis of the partial B X^y decay rate at NNLO 
in renormalization-group improved perturbation theory. A detailed study of the phenomenological 
impact of these effects will be presented elsewhere. 
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